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The genus Dalbergia with c. 250 species has a pantropical distribution. In spite of the high economic and eco- 
logical value of the genus, it has not yet been the focus of a species level phylogenetic study. We utilized ITS 
nuclear sequence data and included 64 Dalbergia species representative of its entire geographic range to pro- 
vide a first phylogenetic framework of the genus to evaluate previous infrageneric classifications based on 
morphological data. The phylogenetic analyses performed suggest that Dalbergia is monophyletic and that 
it probably originated in the New World. Several clades corresponding to sections of these previous classifi- 
cations are revealed. Taking into account that there is not a complete correlation between geography and 
phylogeny, and the estimation that the Dalbergia stem and crown clades are 40.4-43.3 mya and 3.8-12.7 
mya, respectively, it is plausible that several long distance dispersal events underlie the pantropical distribu- 
tion of the genus. 


Systematics 


© 2013 SAAB. Published by Elsevier B.V. All rights reserved. 


1. Introduction 


Understanding evolutionary history of species is a fundamental factor 
for biodiversity assessment and setting conservation priorities (Purvis 
and Hector, 2000), especially when taxa are economically and ecological- 
ly important (Brooks et al., 2006; Margules and Pressey, 2000). Among 
flowering plants, the Leguminosae (Fabaceae) is the third largest family 
with over 19,300 species found in most of the world's vegetation types 
(Lewis et al., 2005). The economic and ecological importance of the le- 
gumes is exceptionally high because they are sources of food, fodder, 
woods, fuels, medicines, and through their ability to enrich soils. 

Because various ecosystems and habitats are under significant 
threat by human activities (Lee and Jetz, 2008; Rockstrom et al., 
2009), monitoring of biodiversity and ecosystem functioning is now 
critical (Yahara et al., 2012). Because of its economic and ecological 
importance, the Leguminosae has been chosen as the focus of a Global 
Diversity Assessment (Global Legume Diversity Assessment; GLDA; 
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Yahara et al., 2013). Within the GLDA, biodiversity assessments of 
multiple legume genera (Bauhinia L., Dalbergia L.f., Desmodium Desv., 
Mucuna Adans., and Vigna Savi.) have been prioritized (Yahara et al., 
2013) implementing Phylogenetic Diversity (PD; Faith, 1992), Function- 
al Diversity (FD; Tilman, 2001) and ecological niche modeling methods. 
One of the genera chosen was Dalbergia, because of its species-richness, 
pantropical distribution, and economic importance. Robust species level 
phylogenies are the starting points for the above-mentioned analy- 
ses, and here we present the first molecular phylogeny of Dalbergia 
(Leguminosae: Dalbergieae). 

The genus Dalbergia is pantropical with around 250 species and 
centers of diversity in Central and South America, Africa, Madagascar 
and Asia (Klitgaard and Lavin, 2005). Several species of Dalbergia 
produce fine timbers of high economic value. These are generally 
known as rosewoods (e.g. D. cochinchinensis Pierre, D. latifolia Roxb., D. 
melanoxylon Guill. & Perr., D. nigra (Vell.) Alleinao ex Benth., D. odorifera 
T.C.Chen and D. sissoo Roxb.,) and are, for example, used in fine carvings. 
More than 50 species of Dalbergia have been demonstrated to fix nitro- 
gen and to have aeschynomenoid type root nodules (Sprent, 2009). 
Several regional revisions, enumerations and flora treatments based 
exclusively on morphological characters exist for Dalbergia (Bentham, 
1860; Carvalho, 1997; Chen et al, 2010; Niyomdham, 2002; Prain, 
1904; Sunarno and Ohashi, 1997; Thothathri, 1987). However, the only 
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attempt at a complete infrageneric classification of the genus is by 
Bentham (1860). 

Combined analyses of a broadly circumscribed monophyletic tribe 
Dalbergieae based on molecular and morphological data (Lavin et al., 
2001), placed Dalbergia in the Dalbergia clade as sister to the genus 
Machaerium Pers. and Aeschynomene L. subgen. Ochopodium. Using 
trnL and nuclear ribosomal DNA (ITS/5.8S) sequence data, Ribeiro et al. 
(2007) resolved Machaerium as more closely related to Aeschynomene 
subgen. Ochopodium Vogel. than to Dalbergia. In that study only 14 spe- 
cies of exclusively New World Dalbergia were included, with the addition 
of the Old World D. sisso from an accession cultivated in Bahia, Brazil 
(Ribeiro et al., 2007). Therefore, a comprehensive phylogenetic study of 
the genus is needed, which includes a representative sample of African, 
Asian and New World species. 

The nuclear ITS region has proven useful in species level phyloge- 
netic studies of various legume genera by many authors (Káss and 
Wink, 1997; Lavin et al., 2000; Wojciechowski et al., 1999). Lavin 
et al. (2001) used ITS sequence data for their broad-scale study of gen- 
era resolved in the Dalbergioid clade and including several Dalbergia 
species. This motivated us to use sequence data of the ITS region for 
this study. We employed a relatively wide taxonomic sampling using 
64 Dalbergia species from its centers of diversity in the Americas, Africa 
and Asia. The objectives of this study are to: (1) provide a first species 
level phylogenetic framework for Dalbergia, (2) evaluate previous 
infrageneric classifications of the genus, and (3) infer the biogeograph- 
ical history of the genus. 


2. Materials and methods 
2.1. Taxon sampling 


Our sampling strategy was to cover the total distribution range of 
Dalbergia and to include representative species from various taxo- 
nomic ranks (subgenus, section and series). In total, 64 species of 
Dalbergia (represented by 66 accessions) were included: 31 species 
from Asia, 21 species from Africa including Madagascar, and 12 spe- 
cies from America (Table 1). Leaf material for DNA extraction was 
obtained from recent field surveys in South East Asia by Tetsukazu 
Yahara and his colleagues, and from E, L, MO and WAG (acronym 
after Thiers, 2013). Twenty-one accessions of Dalbergia species were 
obtained from GenBank. 

In order to evaluate the monophyly of Dalbergia we included 
additional accessions from the Dalbergia clade: Geissaspis cristata 
Wight & Arn., Kotschya africana Endl., Bryaspis humulariodes Gledhill 
and Aeschynomene falcata (Poir.) DC. (Table 1). Designated outgroups 
were selected based on their placement outside of the Dalbergia clade 
in Lavin et al. (2001) and included Adesmia parvifolia Phil., Pterocarpus 
acapulcensis Rose and Arachis correntina (Burkart) Krapov. & W.C.Greg. 
Sequences of all outgroups were retrieved from GenBank. Voucher 
specimens of the newly collected samples in this study were deposited 
in the herbarium of Kyushu University (FU). All newly generated DNA 
sequences in this study are submitted to the DNA Data Bank of Japan 
(DDBJ, AB828599-AB828668). 


2.2. DNA extraction, PCR, and sequencing 


Genomic DNA was extracted from silica gel dried leaves or herbar- 
ium specimens using the DNeasy Plant Mini Kit (Qiagen) following 
the manufacturer's instructions with a modified protocol for herbari- 
um materials. The concentration of genomic DNA was measured with 
a GeneQuant 100 electrophotometer (GE Healthcare, Life Sciences). 
To amplify complete sequence of ITS (ITS1, 5.8S and ITS2) primer 
pairs ITS5-ITS4 (White et al., 1990) were used. In the case of failed 
amplification, for most of the herbarium samples, internal primers 
ITS2 and ITS3 (Baldwin, 1992) were used to amplify the ITS region in 
two pieces with ITS5 and ITSA, respectively. Polymerase chain reactions 


(PCR) were performed in reaction volumes of 10 uL containing 0.2 units 
of ExTaq (TaKaRa) or 0.25 units of MightyAmp DNA Polymerase 
(TaKaRa), and 0.2 mM dNTPs., 10x PCR buffer contains 1.5 mM MgCl, 
0.5-1 uM of each primer pair, and 20 ng genomic DNA. The PCR condi- 
tions were as follows: 2 min for initial denaturation at 95 *C, followed 
by 35 amplification cycles of 45 sec denaturation at 95 *C, 1 min 
annealing at 56 °C, 1 min extension at 72 °C, and a final 10 min exten- 
sion at 72 °C. The PCR products were visualized by performing 0.8% 
agarose gel electrophoresis. Amplified DNA was purified using illustra 
ExoStar Enzymatic PCR and Sequencing Clean-Up kit (GE Healthcare, 
UK) according to the manufacturer's instructions. 

The cycle sequencing reactions were carried out using the ABI 
BigDye Terminator version 3.1 Cycle Sequencing Kit (Applied Biosystems) 
and sequencing reaction products were purified by an ethanol precipi- 
tation method. All DNA sequences were determined in forward and re- 
verse strands using an ABI 3500 DNA sequencer (Applied Biosystems). 
For 10 accessions where direct sequencing yielded electropherograms 
showing polymorphism at several sites, the SSCP (Single strand confor- 
mation polymorphism) method was used (Hayashi, 1991). In this 
method multiple copies of ITS sequences are separated by conformation 
of their molecular structure and weight on an SSCP gel. Distinct single 
bands of SCCP gel were used as PCR templates for sequencing. All for- 
ward and reverse strands for each sequence were assembled, manually 
edited and aligned with MUSCLE (Edgar, 2004). BLAST search (http:// 
blast.ncbi.nlm.nih.gov/) was conducted for all sequences to check for 
possible problems of amplification of contaminant DNA. Given the possi- 
bility of paralogy and/or pseudogenes among ITS sequences, we did not 
include the sequences from SSCP which appeared to be pseudogenes in 
the sequences based upon the presence of indels and high divergence 
of sequences (Bailey et al., 2003). 


2.3. Phylogenetic inferences 


Maximum Parsimony (MP) and Bayesian Inference (BI) were used 
for the phylogenetic analyses. The MP analyses were conducted using 
PAUP* version 4.0b10 (Swofford, 2002). Heuristic searches were 
performed using the tree-bisection-reconnection (TBR) branch swap- 
ping algorithm and all gaps removed from analyses. All characters 
were selected as unweighted and unordered. Branch support was eval- 
uated using the bootstrapping method of Felsenstein (1985) based on 
10000 replicates. 

Bayesian inference was performed using MrBayes version 3.2.1 
(Ronquist et al., 2012). The GIR + I + G model was selected as the 
best-fit model based on the Akaike information criterion (Akaike, 
1974) using the program jModeltest (Darriba et al., 2012; Guindon and 
Gascuel, 2003). Two independent Markov Chain Monte Carlo (MCMC) 
analyses with four simultaneous chains and 8 x 10° generations were 
run. Trees were sampled every 500 generations and the first 4000 
trees were discarded as burn-in. The convergence of MCMC chains 
was visualized with the Tracer program version 1.5 (Rambaut and 
Drummond, 2009) and likelihood scores for the sampled trees were 
inspected. 


3. Results 


The total sequence length of the ITS (ITS1, 5.8S, and ITS2) sequences 
for 64 representatives of Dalbergia and four species from Dalbergia clade 
ranged from 634 to 737 nucleotides sites. When aligned, 362 sites were 
parsimony informative. 


3.1. Phylogenetic inferences 


The maximum parsimony analysis resulted in 4420 most parsimoni- 
ous trees with a length of 1405 steps, consistency index (CI; Kluge and 
Farris, 1969) of 0.429 and retention index (RI; Farris, 1989) of 0.793. 
The MP analysis was not well-resolved, with the strict consensus tree 
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Table 1 
Species list and voucher information for Dalbergia species and four related genera from Dalbergia clade. 
Taxon Author Collector Collection number Continent Country Source 
Dalbergia Lf. 
D. abbreviata Craib Kessler, P.J.A. 2022 Asia Indonesia L 
D. abrahamii Bosser & R. Rabev. Forest 5659 Africa Madagascar L 
D. acariiantha Harms Sally Bidgood & R Abdallah & K Vollesen 1348 Africa Tanzania L 
D. adamii Berhaut Jongkind, C.C.H. 8088 Africa Guinea WAG 
D. afzeliana G.Don Wieringa, J.J. 3844 Africa Gabon WAG 
D. arbutifolia Baker Kayombo, CJ. 1026 Africa Tanzania L 
D. armata E.Mey. Groenendijk, E.M.C. 241 Africa Mozambique WAG 
D. assamica Benth. Ninh, T.N. 11472 Asia Vietnam L 
D. aurea Bosser & R. Rabev. Capuron, R.P.R. 3068 Africa Madagascar WAG 
D. balansae Prain Hu, S.Y. 22805 Asia China L 
D. baronii Baker Rabevohitra, R. 3875 Africa Madagascar WAG 
D. benthamii Prain Hu, S.Y. 21989 Asia China L 
D. bignonae Berhaut Ptsultaf 7983 Africa Benin WAG 
D. bintuluensis Sunarno & Ohashi Sidiyasa, K. 2375 Asia Indonesia L 
D. boehmii Taub. Kayombo, CJ. 6861 Africa Tanzania WAG 
D. bojeri Drake Randrianaivo, R. 694 Africa Madagascar WAG 
D. bracteolata Baker Wohlhauser, S. 60139 Africa Madagascar WAG 
D. candenatensis (Dennst.) Prain Maxwell, J.F. 166 Asia Singapore L 
D. candenatensis (Dennst.) Prain Wightman, G.M. 2701 Asia Australia L 
D. canescens (Elmer) Merr. Soejarto, D.D. 7479 Asia Philippines L 
D. capuronii Bosser & R. Rabev. Capuron, R.P.R. 27807 Africa Madagascar WAG 
D. cearensis Ducke Jardim, J. 1775 America Brazil MO 
D. cochinchinensis Pierre Maxwell, J.F. 1342 Asia Thailand L 
D. cochinchinensis Pierre Nanthavong, K. 643 Asia Laos L 
D. cultrata Pittier Maxwell, J.F. 196 Asia Thailand L 
D. dongnaiensis Pierre Maxwell, J.F. 849 Asia Thailand L 
D. ecastophyllum (L.) Taub. Croat, T.B. 23287 America Belize L 
D. falcata Prain Ambriansyah 448 Asia Indonesia L 
D. floribunda Craib Maxwell, J.F. 159 Asia Thailand L 
D. frutescens var. frutescens (Vell.) Britton Wasum, R. 1585 America Brazil L 
D. frutescens var. tomentosa (Vogel) Benth. Solomon 9409 America Bolivia L 
D. glomerata Hemsl. Sousa, M. 9349 America Mexico L 
D. godefroyi Prain Maxwell, J.F. 65 Asia Thailand L 
D. greveana Baill. Randrianaivo, R. 441 Africa Madagascar L 
D. hancei Benth. Hu, S.Y. 21333 Asia China L 
D. havilandii Prain Ambriansyah 2599 Asia Indonesia L 
D. horrida (Dennst.) Mabb. Reddi, B.V. 25899 Asia India L 
D. hostilis Benth. Jongkind, C.C.H. 7790 Africa Guinea WAG 
D. humbertii R.Vig. Meyers 229 Africa Madagascar L 
D. inundata Benth. Meireles, J. E. 367 America Brazil MO 
D. junghuhnii Benth. Gardette, E. 2177 Asia Malaysia L 
D. kurzii Prain Maxwell, J.F. 286 Asia Thailand L 
D. lactea Vatke De Wilde 6120 Africa Ethiopia WAG 
D. lakhonensis Gagnep. Maxwell, J.F. 444 Asia Thailand L 
D. lanceolaria Lf, Maxwell, J.F. 921 Asia Thailand L 
D. lateriflora Benth. Silva 1057 America Brazil MO 
D. latifolia Roxb. Ambriansyah 963 Asia Indonesia L 
D. maritima RVig. Faliniaina 10 Africa Madagascar L 
D. martinii F.White Brumitt 16910 Africa Zambia L 
D. melanocardium Pittier Wendt, T. 3889 America Mexico L 
D. melanoxylon Guill. & Perr. I FIRIIS 6876 Africa Ethiopia WAG 
D. miscolobium Benth Sano 14516 America Brazil MO 
D. monetaria Lf. Acevedo 7608 America Ecuador U 
D. oliveri Prain Parakosonh, K.P. 161 Asia Laos L 
D. parviflora Roxb. Niyomdham, C. 811 Asia Thailand L 
D. pinnata (Lour.) Prain Huq, A.M. 10430 Asia Bangladesh L 
D. revoluta Ducke s.n. - America Peru MO 
D. rimosa var. foliacea (Benth.) Thoth. FRDU 12 Asia Thailand L 
D. rostrata Hassk. Polak, A.M. 1160 Asia Indonesia L 
D. sandakanensis Sunarno & Ohashi Kessler, P.J.A. 1386 Asia Indonesia L 
D. sericea G.Don Huq, A.M. 10336 Asia Bangladesh L 
D. sissoo DC. Shetty, B.V. 1278 Asia India L 
D. spruceana (Benth.) Benth. Smith 13743 America Bolivia U 
D. stipulacea Roxb. Chin, S.C. 3498 Asia Indonesia L 
D. thorelii Gagnep. Toyama, Tagane & Yahara 1033 Asia Cambodia FU 
D. trichocarpa Baker Service Forestier de Madagascar 21658 Africa Madagascar L 
Aeschynomene falcata (Poir.) DC. Antezana Valera C. 1218 America Bolivia MO 
Bryaspis humularioides Gledhill John Thomas 10097 Africa Liberia MO 
Geissaspis cristata Wight & Arn. Hiep, N.T. 10736 Asia Vietnam MO 
Kotschya africana Endl. Roy E. Gereau 2498 Africa Tanzania MO 
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Fig. 1. The 5075 majority rule consensus tree resulting from Bayesian analysis of the ITS for the members of Dalbergia clade. Numbers along branches are posterior probabilities and 
bootstrap values from Bayesian analysis and Maximum parsimony, respectively. Values <50% were not shown. Bold lines represent clades and subclades with PP >90%. Sequences 


for OTUs with accession number obtained from GenBank. 


showing a basal polytomy. However, major clades are resolved, some of 
which have bootstrap support values greater than 50%. These clades are 
consistent with those found in the Bayesian analyses (Figs. 1 and 2), and 
these have higher posterior probability (PP) than the parsimony boot- 
strap values (though this is usual; Alfaro et al., 2003). The 50% majority 
rule consensus tree of the Bayesian tree revealed better resolution of 
taxa with higher PP support (Figs. 1 and 2). 

Within the Dalbergia clade the West African genus Bryaspis P. A. 
Duvign. is nested within the Aeschynomene clade (Fig. 1). Geissaspis 
cristata and Kotschya africana are also nested within the Aeschynomene 
clade, though their relationships are weakly resolved. 

Dalbergia is monophyletic with two Neotropical species 
(D. miscolobium Benth. and D. spruceana (Benth.) Benth., Clade I) re- 
solved sister to the remaining species, though with relatively weak sup- 
port (0.61 PP). A second New World Clade II consisting of six species 
(Fig. 2) is then sister to all other Dalbergia species (Fig. 2), though again, 
support for this relationship is not exceptionally high (0.88 PP), and nei- 
ther node is resolved in the parsimony strict consensus tree. The species 
in clades I and II belong to sections Dalbergia and Selenolobium Benth. 
sensu Carvalho (1997), and our results show that these sections are 
non-monophyletic. Clade III contains three subclades (II-a, III-b and 
III-c). The well-supported subclade II-a contains D. adamii Berhaut from 
Africa as sister to four members of the Neotropical section Ecastaphyllum 
(P. Browne) Ducke sensu Carvalho (1997). Subclades III-b and Ill-c are 
exclusively Asian. These three subclades are well supported, but relation- 
ships among them are unresolved. Although Clade IV has low PP support 
(PP = 0.54), it brings D. afzeliana G. Don, sampled from Gabon, as sister 
to two subclades (IV-a and IV-b) resolved with PP support of 1 and 0.8, 
respectively. Subclade IV-a includes an African-Madagascan clade (M1) 
as sister to an entirely Asian clade. Subclade IV-b, resolved with relative- 
ly high support (PP = 0.8), reveals a similar pattern to subclade IV-a, 
namely a sister relationship of an African-Madagascan clade with 
a clade consisting entirely of Asian species. Species sampled from 


section Dalbergaria sensu Prain (1904) are in this subclade, though 
with low support (PP = 0.54). Clade V, resolved with high PP sup- 
port (PP = 0.96), contains three major subclades (V-a, V-b and 
V-c) which correspond geographically to the New World (V-a), Africa 
(V-b) and Asia (V-c), respectively, with the exception of D. acariiantha 
Harms from Tanzania which is nested within the Asian subclade V-c 
(Fig. 2). Subclade V-a contains New World species from sect. Triptolemea 
A (Mart. Ex Benth.) Benth. Basally branching to the three well-supported 
subclades in Clade V is a grade consisting of a mixture of New and Old 
World species. 


4. Discussion 
4.1. Phylogenetic relationships within the Dalbergia Clade 


Dalbergioid legumes were recently circumscribed as a pantropical 
group of papilionoids based on molecular and morphological data, 
and initially comprised 44 genera with ca. 1100 species (Lavin et al., 
2001, 2000). Recently, however, three additional genera were placed 
in the Dalbergioid clade: Acosmium Schott s.s. (Cardoso et al., 2012), 
Maraniona C.E. Hughes, G.P. Lewis, Daza & Reynel. (Hughes et al., 
2004), Steinbachiella Harms (Lewis et al., 2012), taking the total to 47. 
Nearly all dalbergioids occur in the Neotropics and Africa. In contrast, 
only about 100 species belonging to six of these 47 genera are distribut- 
ed in tropical Asia. Most of these Asian species belong to Dalbergia 
(Klitgaard and Lavin, 2005; Lavin et al., 2001). In this study we included 
ITS sequence data of four species from Aeschynomene s.l. (Table 1) along 
with Dalbergia clade accessions from GenBank. However, relationships 
among Aeschynomene and segregate genera such as Bryaspis, Geissaspis 
Wight & Arn. and Kotschya Endl. are weakly resolved. Extensive sam- 
pling within Aeschynomene s.s. is required for a better understanding 
of relationships among these genera. 


Fig. 2. The 50% majority rule consensus tree resulting from Bayesian analysis of the ITS for Dalbergia species. Numbers along branches are posterior probabilities and bootstrap 
values from Bayesian analysis and maximum parsimony, respectively. Values <50% were not shown. Bold lines represent clades and subclades with PP >90%. Sequences for 


OTUs with accession number obtained from GenBank. 
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1 1 Dalbergia glomerata Mexico Clade V 
0.53 0.98 Dalbergia melanocardium Mexico V-a 
0.82 Dalbergia frutescens var. tomentosa Bolivia E 1 
Dalbergia brasiliensis EF451076 sect. Tr ipto. 2 
Dalbergia frutescens var. frutescens Brazil Carvalho 19 
a 1 “Dalbergia cearensis Brazil y 
A S l a 0.75 Dalbergia maritinii Zambia 
0.93 |; Dalbergia bojeri Madagascar 
Dalbergia baronii Madagascar 


Af r i ca 0.52] Dalbergia humbertii Madagascar M3 V-b 


0.92 Dalbergia abrahamii Madagascar 
0.99 Dalbergia greveana Madagascar 


A m e r || Cc a 0.97 1 Dalbergia trichocarpa Madagascar 
Dalbergia acariiantha Tanzania 
1 Dalbergia cultrata Thailand 


0.1 Dalbergia sissoo India 
1 1/94 Dalbergia sissoo EF451079 b Wc 
0.97 Dalbergia parviflora Thailand 
0.93 Dalbergia tonkinensis FR854138 
0.96 


Dalbergia rimosa var. foliacea Thailand 
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4.2. Infrageneric classification of Dalbergia 


In 1860, Bentham (1860) divided the 64 species of Dalbergia known 
at that time into six series (Triptolemea Americanae, Triptolemea, Sissoae 
Americanae, Sissoae Gerontogee, Dalbergariae and Selenolobia). In the 
Neotropics, the Brazilian species of Dalbergia were studied intensively 
by Carvalho (1989, 1997). He divided the 41 Brazilian Dalbergia species 
into five sections based on inflorescence and fruit types. In Asia, Prain 
(1904) classified the 86 South East Asian species of Dalbergia into two 
subgenera, five sections and 24 series. Finally, Thothathri (1987) catego- 
rized the 46 Dalbergia species present in the Indian subcontinent into 
four sections and seven series based on androecium and fruit types. In 
general Dalbergia species are morphologically variable, and they possess 
a wide range of habitat preferences which previously made it difficult to 
classify the New and Old World species into natural groups (Bentham, 
1860; Carvalho, 1989; Prain, 1904). Although our ITS analyses only in- 
cluded 64 of the total 250 species in the genus, it contains representa- 
tives from nearly all of the various subgenera, sections and series of 
previous infrageneric classifications. Taking the Bayesian tree (Fig. 2) 
into account, our results revealed that sect. Triptolemea, with cymose 
inflorescences and thin samaroid pods and sect. Ecastaphyllum, with 
racemose or paniculate inflorescences and orbicular to suborbicular or 
reniform fruits, are potentially monophyletic. In contrast, sects. Dalbergia 
and Selenolobium are non-monophyletic (Fig. 2). These results are con- 
gruent with the results of Ribeiro et al. (2007) based on ITS and trnL se- 
quences, and it suggests that inflorescence and fruit types may serve as 
sources of synapomorphies for classifications of Dalbergia as also previ- 
ously noted by Carvalho (1997). Among the Asian Dalbergia species 
only members of sect. Dalbergaria (Prain, 1904) group monophyletically, 
though with low support. Members of this section, with reflexed stan- 
dard petals and stamens usually in two bundles of five, are distributed 
throughout Southeast Asia. 


4.3. Biogeography 


The placement of two Neotropical clades as sister to the remaining 
species in the rest of the genus is suggestive of the New World as the 
area of origin for Dalbergia. However, the low support for these rela- 
tionships (Fig. 2) and the fact that they are not resolved in the parsimo- 
ny strict consensus tree implies that this hypothesis should remain 
preliminary. 

Long distance dispersal (LDD) and vicariance are two processes 
which have been suggested to explain pantropical distributions of vari- 
ous taxa (Bartish et al., 2011). Within the dalbergioid legumes (Lavin 
et al, 2001) the genus Dalbergia has a pantropical distribution with 
high species diversity in Asia (ca. 100 species) as well as in the Neotrop- 
ics and Africa. Using the penalized likelihood method (Sanderson, 2002) 
based on matK sequences, the estimated age of the stem of Dalbergia 
(MRCA of Machaerium copote Dugand and Dalbergia congestifolia Pittier) 
and the crown of Dalbergia (D. sissoo and D. congestifolia) were 40.4-43.0 
mya and 3.8-12.7 mya, respectively, (Matt Lavin, Department of Plant 
Sciences and Plant Pathology, Montana State University, Bozeman, 
Montana; personal communication, see also Lavin et al. (2005, 2004)). 
The Dalbergia stem clade is old because of the age constraint imposed 
on the Machaerium stem clade by a fossil record. The two species of 
crown Dalbergia (D. sissoo and D. congestifolia) in our phylogenetic tree 
are in Clade V (Fig. 2), so the crown age is likely to be an underestimate. 
However it can give a broad time-frame for the age of transoceanic 
crown clades within the Dalbergia clade. Dalbergia species are likely to 
be younger than any of the supposed inter-continental land connections 
or stepping stones (e.g., the North Atlantic; Pennington and Dick, 2004) 
and suggest that multiple migration by LDD across oceans occurred to 
achieve pantropical distribution of Dalbergia. In general Dalbergia species 
have samaroid pods which are adapted to wind dispersal (Van der 
Pijl, 1982). Some species, particularly in section Ecastaphyllum, 
such as D. brownei (Jacq.) Urb., D. monetaria L.f., and D. ecastaphyllum 


(L.) Taub. are adapted to water dispersal (Carvalho, 1997; Croat, 
1978). Additionally, Dalbergia species grow in diverse habitats in- 
cluding moist and dry topical forests, savannas, costal dunes and 
rocky outcrops (Carvalho, 1997; Du Puy and Moat, 2002; Mabberley, 
2008) and contain a variety of woody life forms including shrubs, tree 
and climbing lianas. This diversity of ecology and life forms may have 
contributed to the ability of Dalbergia species to successfully colonize 
and thereby to expand their distribution range. 


4.4. Colonization of Madagascar 


Madagascar is one of the biodiversity hotspots in the world due to 
its high levels of diversity, endemism and its highly threatened biota 
(Yoder and Nowak, 2006). Gondwanan vicariance and/or Cenozoic 
dispersal has been thought to be the two most probable mechanisms 
responsible for the diversity of Madagascar's biota (Queiroz, 2005). Of 
approximately 10,000 native species of higher plants in Madagascar 
about 80% are endemic. Legumes, with ca. 600 endemic species, have 
a major role in Madagascan biodiversity (Puy, 2002). The Malagasy 
Dalbergia species have diverse distribution patterns and grow in differ- 
ent habitats such as evergreen humid forests, coastal forests and rocky 
outcrops (Du Puy and Moat, 2002). Based on the taxonomic revisions 
of Malagasy Dalbergia (Bosser and Rabevohitra, 1996, 2002, 2005), 47 
species are endemic to Madagascar. Despite the diversity of rosewood 
species in Madagascar, they are under high risk of extinction due to 
heavy habitat loss by logging, mining and other human activities 
(Barrett et al., 2010). In this study we included nine endemic species 
(D. abrahamii Bosser & R. Rabev., D. aurea Bosser & R.Rabev., D. baronii 
Baker, D. bojeri Drake, D. capuronii Bosser & R.Rabev., D. greveana Baill., 
D. humbertii R.Vig., D. maritima R.Vig. and D. trichocarpa Baker,) as 
well as D. bracteolata Baker, also native in mainland Africa (Table 1). 
The Malagasy Dalbergia are resolved in three different clades along 
with other species from Africa (Fig. 2). The first is a clade comprising 
D. capuronii and D. maritima as sister to D. boehmii Taub. from Tanzania 
(M1 clade, Fig. 2). The second is D. aurea which is sister to D. lactea 
Vatke from Ethiopia (M2 clade, Fig. 2). The third is a clade of six species 
of Malagasy Dalbergia to which D. martinii F. White from Zambia is the 
sister group with high support (PP — 0.93) (M3 clade, Fig. 2). The 
occurrence of the Malagasy Dalbergia species in three distinct clades 
(M1-3, Fig. 2) suggests at least three independent migrations to the 
island during the Cenozoic. The inference of African species as sister to 
each Madagascan clade suggests that migration may have occurred 
from Africa to Madagascar, but sampling of more African and Malagasy 
species is needed to confirm this. 


5. Conclusions 


Our preliminary molecular systematic analyses of Dalbergia provide 
a first phylogenetic framework for future research and informative re- 
sults regarding previous infrageneric classifications. They also give a 
first biogeographic account of the genus. The results show that the ITS 
region provides sufficient phylogenetic signal among Dalbergia species 
to detect some agreement with classifications based on morphology. 
Additional molecular markers from chloroplast and other nuclear re- 
gions and more exhaustive taxon sampling are needed to be able to elu- 
cidate evolutionary patterns among Dalbergia species and its related 
genera with greater confidence. Moreover, with a more comprehen- 
sively sampled phylogeny, PD methods can be implemented to 
provide assessments of the phylogenetic distinctiveness of Dalbergia 
species that are threatened or at risk of extinction. 
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